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Time-dependent,  thermal  behaviour  of  a  lithium-ion  (Li-ion)  polymer  cell  has  been  modelled  for  electric 
vehicle  (EV)  drive  cycles  with  a  view  to  developing  an  effective  battery  thermal  management  system.  The 
fully  coupled,  three-dimensional  transient  electro-thermal  model  has  been  implemented  based  on 
a  finite  volume  method.  To  support  the  numerical  study,  a  high  energy  density  Li-ion  polymer  pouch  cell 
was  tested  in  a  climatic  chamber  for  electric  load  cycles  consisting  of  various  charge  and  discharge  rates, 
and  a  good  agreement  was  found  between  the  model  predictions  and  the  experimental  data.  The  cell- 
level  thermal  behaviour  under  stressful  conditions  such  as  high  power  draw  and  high  ambient 
temperature  was  predicted  with  the  model.  A  significant  temperature  increase  was  observed  in  the 
stressful  condition,  corresponding  to  a  repeated  acceleration  and  deceleration,  indicating  that  an 
effective  battery  thermal  management  system  would  be  required  to  maintain  the  optimal  cell  perfor¬ 
mance  and  also  to  achieve  a  full  battery  lifesapn. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  (Li-ion)  polymer  batteries  potentially  have  much 
higher  energy  density  ( ~  400  Whl-1)  compared  to  other  battery 
technologies,  and  they  are  already  replacing  nickel-metal  hydride 
(NiMH)  batteries  for  portable  devices  such  as  laptop  computers  and 
smart  phones,  and  also  for  hand-held  devices  in  the  power  tool 
industry  [1].  The  high  energy/power  density  of  Li-ion  batteries  has 
also  made  them  an  attractive  power  source  in  electric  vehicles  (EV) 
and  hybrid  electric  vehicle  (HEV),  as  they  can  provide  longer 
driving  range  and  higher  acceleration.  A  typical  battery  pack  for  the 
EV  application  is  composed  of  a  large  number  of  cells  to  provide  the 
required  power,  and  the  compactness  of  battery  packs  poses 
a  challenge  for  effective  thermal  management  [1,2]. 

The  major  obstacles  to  designing  a  scaled-up  of  Li-ion  batteries 
for  EV  applications  have  been  safety,  performance  and  overall  cell 
life.  The  safety  issue  concerns  the  thermal  runaway  of  battery  cell, 
and  the  prediction  and  prevention  of  the  propagation  of  thermal 
runaway  have  been  studied  by  many  researchers.  A  Li-ion  cell 
generates  heat  during  both  charge  and  discharge,  and  the  heat 
generation  in  the  battery  can  increase  sharply  leading  to  over¬ 
heating  under  stressful  conditions  such  as  high  power  draw  and 
high  ambient  temperature  while  overheating  can  also  be  caused  by 
defects  in  individual  cells.  If  the  cells  are  not  thermally  managed, 
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the  thermal  runaway  propagates  into  an  entire  battery  pack. 
Various  battery  thermal  management  systems  have  been  devel¬ 
oped  for  air-  and  liquid-cooling  to  prevent  the  thermal  runaway 
propagation  without  over-designing  the  cooling  system  [3]. 

The  battery  cell  temperature  also  affects  the  performance,  reli¬ 
ability  and  lifespan.  High  battery  cell  temperature  can  cause 
a  shortening  of  the  lifespan  of  the  Li-ion  battery  due  to  increased 
degradation  of  the  battery  cell.  It  was  found  that  every  degree  of 
cell  temperature  rise  reduces  the  lifespan  of  the  Li-ion  battery  by 
approximately  two  months  in  an  operating  temperature  range  of 
30-40  °C  [4].  To  achieve  a  full  lifespan,  the  cell  temperature 
difference  should  be  kept  below  5  °C  [5,6]  as  large  temperature 
non-uniformity  in  the  battery  pack  affects  adversely  overall  cell 
lifespan  as  well  as  available  cell  capacity  [7-10].  Therefore, 
a  battery  thermal  management  system  enabling  effective  temper¬ 
ature  control  becomes  essential  for  safety,  performance,  reliability 
and  overall  cell  life  [2,11,12].  For  better  design  of  thermal 
management  system,  understanding  the  real  time  behaviour  of  Li- 
ion  cells  in  vehicle  operating  conditions  is  critical,  and  modelling  its 
thermal  behaviour  in  various  EV  drive  cycles  can  help  to  design 
such  battery  thermal  management  systems  [1]. 

Most  of  the  previous  studies  have  modelled  the  Li-ion  cell 
thermal  behaviour  separately  for  charge  or  discharge,  and  have  not 
considered  the  combined  effect  of  them,  typically  present  in  vehicle 
load  conditions  [13].  Accurate  modelling  of  the  electro-thermal 
behaviour  of  Li-ion  cells  for  dynamic  electric  loading  conditions  is 
a  key  requirement  to  improve  battery  thermal  management 
system.  In  this  study,  a  three-dimensional,  fully  coupled  electro- 
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Cover  for  controlling  forced  convection  around  cells 


Fig.  1.  a)  Schematic  of  Li-ion  cell  test  setup,  b)  instrumented  cell  in  climatic  chamber,  and  c)  location  of  thermocouples  in  mm. 


thermal  cell  model,  which  is  capable  of  modelling  vehicle  load 
conditions,  is  developed.  A  high  energy  density  Li-ion  polymer 
pouch  cell  is  tested  at  various  charge  and  discharge  rates  to  provide 
model  characterisation  data  and  also  other  relevant  electrical  and 
thermal  data  for  validation.  The  current  work  focuses  on  electro¬ 
thermal  modelling  and  its  use  in  battery  thermal  management 
studies  for  EV/HEV  applications.  This  model  allows  for  scalability 
with  accuracy  and  reasonable  simulation  time  which  is  of  prime 
importance  while  modelling  full  vehicle  battery  models. 

2.  Experimental  setup 

In  this  study,  a  high  energy  density  Li-ion  polymer  pouch  cell 
was  tested  to  assess  the  electro-thermal  cell  model.  The  Li-ion  cell 
has  a  graphite  anode  and  a  Li[NiMnCo]02  cathode  with  a  nominal 
capacity  of  20  A  h.  The  schematic  of  the  experimental  setup  is 

Table  1 

Cell  test  procedure.  rtest  is  the  test  temperature. 


Initialization  step 


Temperature  (°C) 

Procedure 

Value 

Limit  criteria 

25 

Discharge 

1C 

l/min  <  3.0V 

25 

Break 

10  min 

25 

Charge 

1C 

Um ax  >  4.1V 

/min  <  1A 

Measurement  step 

Test 

Break 

>  40  min 

Tell  =  Test 

Test 

Discharge 

various  C-rates 

//min  <  3.0V 

25 

Break 

30  min 

25 

Charge 

IC 

//max  >  4.1V 

/min  ^  1A 

shown  in  Fig.  1.  A  climatic  chamber  was  used  to  provide  controlled 
thermal  conditions  during  the  test  (Fig.  1  a).  A  box  was  installed 
around  the  cells  to  ensure  well-defined  free  convective  cooling  and 
to  protect  the  cells  from  forced  convection,  which  would  adversely 
affect  the  results.  Fig.  1  b  shows  the  Li-ion  cells  inside  the  box.  The 
cell  was  prepared  with  10  thermocouples  to  measure  the  temper¬ 
ature  evolution  on  the  cell  surface  during  the  test  run.  The  locations 
of  the  thermocouples  are  shown  in  Fig.  1  c. 

The  Li-ion  cell  was  tested  using  one  of  the  Tata  Motors  European 
Technical  Centre  (TMETC)  cell  testing  program  within  operating 
limits  of  the  cell.  The  test  cycle  consisted  of  several  constant  current 
charge  and  discharge.  Temperature  was  kept  constant  during  the 
measurement  at  a  reference  temperature  (Tref).  The  electrical  and 
thermal  measurements  at  different  reference  temperatures  and 
discharge  rates  were  recorded.  The  cell  was  discharged  at 
a  constant  discharge  rate  until  the  lower  cut-off  limit  of  3  V  was 
reached.  Each  discharge  rate  measurement  required  initialisation 
steps  to  precondition  the  cell  before  the  measurement.  These  steps 
were  to  ensure  that  the  prescribed  electrical  and  thermal  initial 


Fig.  2.  Charge/discharge  test  cycle. 
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Table  2 


Material  properties  used  in  the  model,  p  is  density,  k  is  thermal  conductivity,  ke  is 
electrical  conductivity,  and  cp  is  specific  heat. 


P  [kg  m  3] 

k  [W  mk1] 

ke  [S  mm  l] 

cP  [J  kg  1  K  M 

Pouch 

1.15 

0.16 

0.1 

1900 

Positive  Terminal 

2.7 

238 

38300 

900 

Negative  Terminal 

8.36 

398 

63300 

385 

condition  for  each  measurement  was  achieved.  Table  1  shows  the 
procedure  for  initialisation  and  measurement  requirements  in 
detail. 

Fig.  2  shows  the  electric  loading  cycle.  The  section  AB  corre¬ 
sponds  to  a  discharge  current  of  5C.  The  constant  current  discharge 
was  followed  by  a  10  min  break  with  no  electric  load.  The 
remaining  initialisation  step  before  the  next  constant  discharge 
consists  of  a  1C  constant  current  charge  (section  EF)  and  a  constant 
voltage  section  (section  FG)  where  the  current  gradually  reduces 
below  1  A.  A  40  min  break  is  then  given  for  the  relaxation  processes 
in  the  cell  chemistry  before  the  next  discharge  rate  measurement. 
This  process  of  initialisation  was  performed  between  the  discharge 
measurements:  sections  HI,  MN  and  RS  correspond  to  0.5C,  1C,  and 
3C  constant  current  discharge  measurements,  respectively.  This 
whole  process  is  repeated  for  different  reference  temperatures 
(Tref),  and  thus  providing  the  cell  characterisation  for  various 
discharge  rates  and  reference  temperatures.  The  cell  temperature 
and  voltage  during  charge  and  discharge  in  the  entire  test  cycle 
were  recorded  and  compared  with  the  the  electro-thermal  model 
prediction. 


modelled  using  a  simplified  polatisation  characteristics  of  the 
electrodes,  and  Poisson  equations  for  the  potentials  of  the  positive 
and  negative  electrodes  (Vp  and  Vn)  are  derived  in  conjunction  with 
Ohm’s  law.  The  cell  terminal  voltage  can  then  be  obtained  as 
Vceii  =  Vp  -  Vn.  A  modified  form  of  the  polarisation  expression 
proposed  by  Newman  and  Tiedemann  [15]  is  used  in  this  study  to 
estimate  the  current  density  (J).  An  additional  temperature 
dependency  term  is  included  in  the  polarisation  expression  to 
account  for  the  Seebeck  effect.  In  the  following  equation,  the 
current  density  (J)  is  expressed  as  a  function  the  cell  terminal 
voltage  (Vceii)  and  temperature: 

J  =Y(vcen-U-b(Ti-Tref)),  (1) 

where  Y  is  the  conductance,  and  U  is  the  equilibrium  voltage.  Vand 
U  are  estimated  using  the  measurement  data  of  Vceii  and  J.  In  the 
test,  Vceii  is  measured  as  a  function  of  depth-of-discharge  (DOD), 
electric  current  and  temperature.  The  first  term  on  the  right  hand 
side  is  the  Ohm’s  law  term,  and  the  current  density  is  defined  as 
J  =  I/A,  where  /  is  the  current  normal  to  the  active  material,  and  A  is 
the  active  material  area.  The  last  term  accounts  for  the  Seebeck 
effect,  where  T\  is  the  cell  temperature,  and  Tref  is  the  ambient 
temperature.  The  Seebeck  coefficient  b  is  to  be  determined  from  the 
test  data. 

The  local  heat  generation  rate  of  the  Li-ion  cell  is  then  estimated 
during  charge/discharge  cycle  based  on  the  various  losses  in  the  cell 
[16].  The  heat  generation  equation  developed  by  Bernardi  etal  [17] 
is  used  in  this  study: 


3.  Coupled  electro-thermal  model 

The  Li-ion  cell  used  in  this  study  is  comprised  of  a  pouch  case, 
electrode  plates,  electrolyte,  and  separators.  To  model  the  thermal 
behaviour  of  a  battery  cell  accurately,  the  configuration,  geometry, 
physical,  chemical  and  electrochemical  properties  need  to  be 
known.  In  our  electro-thermal  model,  the  cell  is  characterised  by 
a  structure  formed  with  a  homogeneous  layer,  representing  the 
anode,  cathode  and  separator.  The  local  heat  generation  in  the 
battery  cell  due  to  the  electrochemical  reactions  and  the  mass 
transfer  of  ions  in  the  electrolyte  are  characterised  by  local  internal 
resistance  and  the  current  densities.  The  performance  of  a  cell  is 
influenced  by  the  charge/discharge  rate,  the  ohmic  resistance,  the 
activation  polarisation  or  charge  transfer  at  the  electrodes,  the 
concentration  polarisation  in  the  electrolyte  near  the  electrodes, 
and  the  cell  temperature. 

In  this  study,  a  three-dimensional  mathematical  model  is  used 
to  predict  the  electrical  and  thermal  performance  of  the  Li-ion  cell. 
The  potential  and  current  density  distribution  on  the  positive  and 
negative  electrodes  are  calculated  using  a  PDE  method  proposed  by 
Kwon  et  al  [14].  The  charge  conservation  on  the  electrodes  are 


Q  =  l(u-VceU-T^j,  (2) 

where  I(U  -  Vcen)  estimates  irreversible  heat  component  due  to 
deviation  of  cell  voltage  from  the  equilibrium  voltage,  and  I(T[dU/ 
dT )  accounts  for  the  reversible  heat  generation  due  to  entropy 
change.  dU/dT  expresses  the  temperature  dependence  of  the  equi¬ 
librium  voltage.  The  heat  generation  due  to  the  ohmic  resistance  in 
the  cell  is  also  considered  in  the  study  [18].  Please  note  that  the  cell 
terminal  voltage  (Vcen  =  Vp  -  Vn)  is  obtained  from  the  Poisson 
equations,  as  explained  earlier. 

In  order  to  model  the  thermal  behaviour  of  the  Li-ion  cell,  the 
electrical  model  is  combined  with  an  energy  balance.  The  unsteady 
three-dimensional  heat  conduction  equation  can  be  written  as 
[19-21]: 


9  T 

pCpdF 


0  (  0T\  0  (  0T\ 

aA'W  +0^V y9yJ 


+  Q? 


(3) 


where  p  is  the  density,  cp  the  heat  capacity,  and  k  the  effective 
thermal  conductivity,  and  Q  the  heat  generation  in  the  Li-ion  cell 
(Eq.  (2)). 


Fig.  3.  Li-ion  cell  characteristics  at  several  ambient  temperatures,  a)  Capacity,  and  b)  resistance  variation  with  Tre f. 
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Fig.  4.  Comparison  of  electro-thermal  model  prediction  with  experimental  test  results  at  Tre f  =  25  °C.  a)  Average  voltage,  and  b)  average  temperature  versus  DOD  at  various 
discharge  rates. 


A  convection  thermal  boundary  condition  is  applied  [17]  at  the 
cell  surface,  to  consider  the  heat  transfer  (Qc)  with  the 
environment, 

Oc  =  h(r-rref),  (4) 

where  h  =  25  W  m_2I<1  is  the  convective  heat  transfer  coefficient. 
The  reference  temperature  is  set  to  be  Tre f  =  25  °C  to  reflect  the 
climatic  chamber  test  temperature. 

Eqs.  (2)  and  (3)  can  predict  the  electro-thermal  behaviour  at 
every  point  in  the  Li-ion  cell.  The  heat  release  from  the  cell,  and 
subsequently  the  cell  temperature  are  solved  in  both  space  and 
time.  A  three-dimensional  finite  volume  method  is  used  to  solve 
the  electro-thermal  model,  and  the  cell  volume  is  discretised  as 
a  cartesian  mesh  in  all  three  directions.  In  the  finite  volume  model, 
the  electrodes,  separators  and  electrolyte  are  modelled  at  a  macro 
level,  and  the  volume-averaged  material  properties  of  all  the  cell 
components  are  used  [21  ].  The  enveloping  pouch  is  also  considered 
to  model  a  realistic  heat  transfer  to  the  environment.  The  material 
properties  of  the  pouch  and  the  tab  used  in  the  simulation  are  listed 
in  Table  2.  The  density,  thermal  conductivity  and  heat  capacity  of 
the  cell  components  are  assumed  to  be  uniform  throughout  the 
battery  and  to  remain  constant  within  a  known  range  of 
temperature. 

Y  and  U  in  Eqs.  (1)  and  (2)  are  dependent  on  the  depth-of- 
discharge  (DOD)  of  the  cell.  Y  and  U  values  are  estimated  at 
various  DOD  using  measurement  data  of  cell  voltage  versus  current 
density,  where  Y  is  the  inverse  of  the  slope,  and  U  is  the  intercept. 
Then,  their  DOD  dependency  is  modelled  following  Gu  [22].  The 
conductance  is  given  by: 


y(DOD)  =  ^ak(DOD)kexp 

k 


Eact 

~R~ 


(5) 


where  R  is  the  universal  gas  constant  (8.314  J  mol-1  K-1)).  Arrhe¬ 
nius  expression  [23]  is  used  to  consider  the  temperature  depen¬ 
dence  of  the  electrochemical  model  property.  The  activation 
energy,  Eact,  is  to  be  estimated  from  the  test  data.  Similarly,  the 
equilibrium  voltage  is  given  by 


If  (DOD)  =  ]TMDOD)m, 

m 


(6) 


where  a k  and  bm  are  the  fitting  coefficients  for  Yand  U,  respectively. 
In  this  study,  all  fitting  coefficients  could  represent  up  to  8th-order 
polynomials,  and  these  coefficients  are  determined  from  discharge 
test  data  for  various  temperatures  and  C-rates.  Similarly,  dU/dT  in 
Eq.  (2)  is  also  characterised  with  8th-order  polynomials. 

Now,  the  values  for  b  and  Eact  can  be  determined  using  discharge 
test  data.  From  Eqs.  (1),  (5)  and  (6)  we  have 


Fig.  5.  Comparison  of  electro-thermal  model  prediction  with  experimental  test  results, 
a)  Current,  b)  SOC,  c)  average  voltage,  and  d)  temperature  at  P01. 
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Y(DOD,  Tj) 


5>k,i(  DOD)k 


y^qk(DOD)kexp 

k 


Fact 

~R~ 


different  ambient  temperatures  (Ire f).  Fig.  3  shows  that  the  usable 
cell  capacity  increases  significantly  with  Fref  while  the  internal 
resistance  of  the  Li-ion  cell  decreases.  The  cell  has  the  maximum 
usable  capacity  in  a  range  of  25  <  Fref  <  45  °C,  where  the  internal 
resistance  remains  small. 


(J(DOD)  +  b(Tj  -  Tref)  =  ^Bmi(DOD)m 

m 

=  £MDOD)m+b(Ti-rref),  (8) 

m 

where  Akij  and  Bm,i  are  fitting  coefficients  for  a  given  temperature  T\. 
From  Eqs.  (7)  and  (8),  £act  and  b  can  be  deduced 


fEiA,i(DOD)kA 

V  EkQk(DOD)k  J 


-act 


1  /  1 


E  I  Fpef 


(9) 


where  £act  is  the  slope  of  the  linear  fitting  of 
ln(EkAc,i(DOD)k/Ekak(DOD)k)  with  respect  to  1/R(l/Tref-1/Ti). 
and  b  is  the  slope  of  the  linear  fitting  of  (Bo  -  bo)  to  (Tj  -  Tref).  Test 
data  at  various  temperatures,  charge  and  discharge  rates  are  used 
for  cell  parameter  characterisation. 


4.  Results  and  discussion 


Before  describing  the  details  of  the  electro-thermal  model 
validation  and  prediction,  key  measured  cell  characteristics  are 
summarised  below.  The  change  in  cell  capacity  and  resistance 
during  discharge  with  a  C-rate  of  0.5C  was  measured  at  five 


4.2.  Electro-thermal  model  validation 

A  three-dimensional  transient  electro-thermal  model  is  used  to 
predict  unsteady  temperature  distribution  on  the  Li-ion  cell  for 
a  given  electric  load  cycle  and  thermal  boundary  condition.  The 
proposed  electro-thermal  model  is  correlated  with  test  data  for 
various  discharge  rates,  and  the  DOD  dependency  of  the  electrical 
properties  in  Eqs.  (5)  and  (6)  are  determined. 

The  electro-thermal  model  predictions  are  compared  with 
experimental  test  results  for  several  discharge  rates.  Fig.  4  shows 
the  average  cell  voltage  and  temperature  variation  during  discharge 
at  Tref  =  25  °C.  The  electrical  behaviour  prediction  compares  very 
well  with  the  experimental  test  results  over  a  wide  range  of  DOD. 
The  sudden  voltage  drop  observed  at  the  end  of  discharge  (at  high 
DOD)  is  correctly  predicted  in  the  model.  The  average  cell 
temperature  is  also  in  good  agreement  with  the  test  results.  It  can 
be  seen  that  there  is  a  steady  increase  in  temperature  during  the 
constant  current  discharge.  As  the  electrical  load  increases,  the  heat 
generation  (given  in  Eq.  (2))  and  hence  the  cell  temperature 
increase.  As  expected,  the  largest  temperature  rise  occurs  for  the 
highest  discharge  rate  (higher  electrical  load)  due  to  the  increased 
heat  generation  from  the  cell. 

A  Li-ion  cell  test  cycle  with  a  combination  of  various  charge  and 
discharge  rates  is  then  used  for  validation  of  charge  and  discharge 


Fig.  6.  Electro-thermal  model  temperature  contour  on  cell  face.  A  circle  shows  the  temperature  variation  on  the  side  wall,  a)  Point  U:  5C  discharge,  b)  Point  V:  1C  charge,  c)  Point  W: 
constant  voltage  step,  d)  Point  X:  0.5C  discharge,  e)  Point  Y:  1C  discharge,  and  f)  Point  Z:  3C  discharge. 
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behaviour  modelling.  Fig.  5  shows  a  comparison  between  the 
experimental  test  data  and  the  electro-thermal  model  prediction. 
Various  cell  properties  are  gathered  during  the  test,  and  the  time 
variation  of  current,  SOC,  average  cell  voltage,  and  temperature  are 
shown  in  the  figure.  A  measuring  point  (P01)  is  located  near  the 
negative  tab  (see  lc).  The  overall  electrical  and  thermal  results  from 
both  test  and  electro-thermal  model  are  in  good  agreement.  The 
initial  cell  temperature  is  the  same  as  the  ambient  temperature  of 
Tref  =  25  °C.  At  the  beginning  of  the  cycle,  the  maximum  applied 
discharge  current  of  5C  causes  a  sharp  increase  in  the  cell 
temperature  by  more  than  25  °C  within  approximately  7  min,  and 
the  amount  of  this  temperature  increase  is  correctly  predicted  in 
the  model.  The  immediate  voltage  drop  during  the  5C  discharge  is 
also  correctly  captured.  The  voltage  drop  is  followed  by  a  relaxation 
phase  with  a  sudden  voltage  rise.  This  is  because  the  internal 
resistance  decreases  with  the  temperature  rise  due  to  faster 
kinetics,  hence  resulting  in  an  increased  cell  voltage.  Following  the 
residual  temperature  build  up  from  the  5C  discharge  and  relaxation 
phase,  the  temperature  remains  below  30  °C  during  the  charge.  A 
similar  trend  of  voltage  and  temperature  response  is  also  seen  for 
an  applied  discharge  current  of  3C  at  the  end  of  the  cycle,  and  the 
instantaneous  increase  in  temperature  is  about  15  °C.  A  modest 
temperature  increase  during  1C  and  0.5C  discharge  is  also  well 
predicted  in  the  model,  displaying  a  much  lower  slope  of  the 
temperature  history.  As  expected,  SOC  decreases  linearly  in  time 
during  constant  current  charge  and  decreases  during  discharge. 

Fig.  6  shows  the  model  prediction  of  temperature  contours  on 
the  cell  surface  during  the  test  cycle.  Several  phases  in  the  test  cycle 
are  chosen  to  display  the  battery  thermal  characteristics  during 
both  charge  and  discharge.  U,  V,  W,  X,  Y,  and  Z  points  during  the  test 
cycle  are  depicted  in  Fig.  2.  Temperature  on  the  positive  terminal  is 
higher  than  the  negative  terminal  due  to  the  lower  electrical 
conductivity  of  the  positive  electrode  as  shown  in  Table  2.  For 
higher  discharge  rates  of  3C  and  5C,  the  current  flow  is  significantly 
high  in  the  tab  area  due  to  an  increased  current  density  (in  Eq.  (1)), 
hence  high  temperature,  particularly  near  the  positive  tab.  The 
temperature  variation  on  the  cell  progressively  reduces  for  lower 
discharge  rate  as  shown  in  the  contour  plot.  The  three-dimensional 
model  is  able  to  predict  the  temperature  variation  on  the  cell 
surface.  A  similar  temperature  variation  on  the  tabs  and  cell  were 
observed  in  the  recent  cell  test  thermal  images  of  Kim  et  al  [18]. 

The  results  reported  in  this  Section  has  provided  a  good  vali¬ 
dation  for  the  current  electro-thermal  model.  An  excellent  agree¬ 
ment  is  found  between  the  Li-ion  cell  model  prediction  and  the  test 
data  for  charge  and  discharge,  which  clearly  demonstrates  the 
capability  of  the  current  Li-ion  cell  model  to  predict  electro¬ 
thermal  behaviour  of  the  cell  for  different  loading  cycles. 

4.2.  Electro-thermal  model  prediction  for  vehicle  operating 
conditions 

The  electro-thermal  model  is  now  applied  to  predict  the  Li-ion 
cell  behaviour  for  vehicle  operating  conditions.  Two  vehicle  oper¬ 
ating  conditions  for  the  electro-thermal  simulation  are  considered 
in  this  section.  The  first  one  (Cycle  A)  is  the  drive  cycles  during  a  hot 
environment  test  of  an  actual  electric  vehicle.  From  the  TATA  EV 
(prototype  without  cooling)  hot  environment  tests  at  an  proving 
ground  during  the  summer,  a  gradual  temperature  increase  was 
reported  for  12%  grade  climb  at  40  km  h-1  [24].  It  was  also  observed 
that  at  the  end  of  the  battery  capacity  there  was  a  significant 
increase  in  battery  temperature,  indicating  changes  in  internal 
resistance.  The  second  one  (Cycle  B)  is  for  a  more  severe  test,  where 
the  vehicle  was  repeatedly  accelerated  from  20  km  h-1  to 
70  km  h-1  and  then  braked  back  to  20  km  h_1.  In  this  case,  the 
battery  temperature  traces  exhibited  a  small  ‘saw-tooth’  effect  of 


around  0.2-0.5  °C  due  to  the  transient  peaks  in  the  battery  charge 
and  discharge  [24]. 

Fig.  7  shows  the  cell-level  electro-thermal  prediction  for  Cycle  A. 
The  simulation  is  performed  for  a  constant  discharge  cycle  of  1.2C 
(24  A),  representing  an  electrical  loading  for  a  40  km  hr1,  12% 
(ascending  slope)  climb  vehicle  operating  condition.  Fig.  7  shows 
the  time  variation  of  current,  SOC,  average  cell  voltage,  and 
temperature  on  the  cell  surface  for  Cycle  A.  Three  measuring  points 
(P01,  P03,  and  P10)  are  located  near  the  negative  tab  (P01),  and 
positive  tab  (P03  and  P10).  The  SOC  prediction  shows  that  the  cell 
would  completely  discharge  in  about  3000  s  as  would  be  expected 
for  a  cell  of  20  A  h  nominal  capacity.  The  temperature  increases  at 
the  measuring  points  are  gradual,  and  similar  to  the  one  observed 
in  the  vehicle  level  hot  environment  tests.  The  temperature  rise 
time  constant  is  approximately  300  s.  There  is  an  additional 
increase  in  the  temperature  at  the  end  of  discharge,  which  corre¬ 
sponds  to  an  increased  resistance. 
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Fig.  7.  Electro-thermal  model  prediction  for  a  mild  electric  vehicle  operating  condition 
(Cycle  A),  a)  Current,  b)  SOC,  c)  average  voltage,  and  d)  temperature  at  P01,  P03 
and  P10. 
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The  cell-level  electro-thermal  model  is  then  applied  to  a  more 
severe  vehicle  operating  condition  (Cycle  B),  where  a  vehicle  is 
repeatedly  accelerated  and  decelerated  between  20  km  h-1  and 
70  km  h_1.  Each  cycle  of  30  s  is  repeated  for  a  total  of  150  s  to 
understand  the  accumulated  thermal  behaviour.  The  discharge  step 
that  represents  acceleration  was  simulated  for  1.7C  ( -34  A),  and  the 
charge  step  that  represents  braking  was  simulated  for  0.83C  (17  A). 
Fig.  8  shows  the  time  variation  of  current,  SOC,  average  cell  voltage, 
and  temperature  at  three  measuring  points  (P01,  P03,  and  P10)  on 
the  cell  surface  for  Cycle  B.  A  small  saw-tooth  effect  is  clearly  seen 
in  the  cell  temperature  response  at  the  charge/discharge  transients, 
and  a  similar  feature  was  observed  on  TATA  EV  (prototype  without 
cooling)  during  the  vehicle  level  hot  environment  tests  [24].  It  is 
found  that  the  temperature  increases  by  more  than  1  °C  during 
a  aggressive  cycle  of  150  s  at  a  modest  ambient  temperature.  This 
would  mean  a  substantial  temperature  increase  for  a  standard  1  h 


Fig.  8.  Electro-thermal  model  prediction  for  a  severe  electric  vehicle  operating 
condition  (Cycle  B).  a)  Current,  b)  SOC,  c)  average  voltage,  and  d)  temperatures  at  P01, 
P03  and  P10. 


aggressive  hot  environment  test.  Considering  the  fact  that  optimum 
Li-ion  cell  performance  is  found  at  around  25  °C  cell  operating 
point  (see  Fig.  3),  this  would  mean  that  for  such  aggressive  drive 
cycles  the  cells  would  potentially  need  an  optimum  active 
temperature  control  [2]. 

5.  Conclusions 

A  three-dimensional  electro-thermal  cell  model  is  developed  for 
predicting  the  electrical  and  thermal  behaviour  of  a  Li-ion  polymer 
cell.  The  fully  coupled  model  has  been  implemented  based  on 
a  finite  volume  method.  The  DOD  dependency  of  the  cell  is  char¬ 
acterised  using  the  experimental  test  data.  The  Li-ion  cell  test  data 
were  utilised  to  validate  the  electro-thermal  model.  The  electro¬ 
thermal  model  predicts  the  cell  temperature  and  voltage  magni¬ 
tudes  accurately  for  the  test  load  cycle.  The  electro-thermal  model 
is  then  used  to  predict  the  cell-level  performance  for  two  different 
vehicle  drive  cycles:  one  for  a  mild  EV  operating  condition,  and  the 
other  for  a  more  severe  operating  condition.  It  is  found  that 
a  significant  temperature  rise  occurs  during  the  severe  drive  cycle, 
and  an  effective  temperature  control  is  required  to  maintain  the 
optimal  cell  performance.  The  model  has  shown  the  potential  to  be 
used  in  battery  thermal  management  studies  for  EV/HEV  applica¬ 
tions.  The  electro-thermal  data  and  insights  on  cell  surface 
temperature  variation  would  aid  better  design  of  battery  thermal 
management  systems.  In  a  separate  study,  the  proposed  electro¬ 
thermal  cell  model  is  being  implemented  to  predict  cell  thermal 
performance  estimations  in  battery  management  system  and 
vehicle  drive  cycle  modelling  applications.  The  current  work  is  also 
being  extended  to  model  the  cell  electrical  dynamics  and  their 
effect  on  the  voltage  and  temperature. 
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